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Abstract 



Dilepton production associated with minijets is calculated in ultrarela- 
tivistic heavy ion collisions using the first order approximation of the dilepton 
fragmentation functions of quarks and gluons. The full QCD evolution of the 
fragmentation functions is also studied. We find that the dilepton pairs from 
the fragmentation of minijets are comparable to direct Drell-Yan at y/s = 200 
AGeV for small dilepton invariant mass M ~ 1-2 GeV/c 2 while dominant 
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I. INTRODUCTION 



In the search for a quark gluon plasma (QGP) in ultrarelativistic heavy ion collisions, 
electromagnetic signals are considered good probes of the dense matter JIJ. Because of 
the large electromagnetic mean free path, leptons and photons produced by interacting 
(anti) quarks inside QGP can easily escape the hot and dense matter and carry the informa- 
tion of the system to the detector. Recent developments in parton transport phenomenology 
indicate that dileptons and photons could also reveal the dynamics of the early evolution 
of the a dense parton system However, like for all the proposed QGP signals, back- 

ground must be understood and subtracted in order to distinguish the true features of a 
QGP. In general, there are two kinds of background sources for the thermal electromagnetic 
signals. One comes from the evolution of the hadronic phase and the decays of the produced 
hadrons. The other one is due to the initial parton scatterings at the very earliest stage 
of the heavy ion collisions. For dilepton production, the latter one is usually referred to as 
Drell-Yan (DY) processes ||. 

In the lowest order, 0(a 2 ), the DY processes are simply quark-antiquark annihilations. 
First order contributions 0(a 2 a s ) in perturbative QCD (pQCD), originating from initial 
state radiation and virtual corrections, give rise to about the same amount of dilepton pro- 
duction as in the lowest order, which is often characterized by a so-called "iT-factor" of about 
2 These corrections are also responsible for large pt tails of the dilepton transverse 

momentum spectrum. At small p?, summation over the initial state soft gluon radiations 
generates a Sudakov form factor regularizing the perturbative low px production f8fl. By 
now, there also exist matrix element calculations of the second order pQCD contributions, 
C{a 2 a 2 ), to the ^-factor [|. 

In this paper, we will investigate dilepton production associated with minijet final state 
radiation in heavy ion collisions at collider energies. It is expected that at energies \fs > 200 
AGeV, minijets [(anti)quarks and gluons with pt ~ few GeV/c] are produced abundantly 
via multiple semihard scatterings. These minijets have important contributions to particle 
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production, transverse energy and overall evolution of the formed quark-gluon system 



fl"3|| . Therefore, it would be interesting to study dilepton bremsstrahlung from the initially 
produced minijets. Especially, it is important to know whether the dileptons associated 
with minijets could compete with the lower order DY processes at midrapidity and small 
invariant dilepton masses M ~ 1 — 3 GeV/c 2 , where a window for observing the thermal 
dileptons is expected 

Rather than strictly applying the almost complete 0{a 2 a 2 s ) results from the matrix 
element calculations |§ , we take a different and simpler approach by calculating the dilepton 
fragmentation functions of the final state minijets. Unlike in the real photon fragmentation 



functions [0,0 , the relatively large invariant masses M ^> A of the dileptons fix the lower 
limit of the momentum scale of the QCD radiation processes. This makes the problem 
calculable in pQCD. In the leading logarithm approximation and in an axial gauge ||16|| , 



the dilepton fragmentation functions can be calculated up to all orders in pQCD. Using the 
obtained fragmentation functions to convolute with minijet cross sections, we then compute 
the contribution to the dilepton production from the final state radiation of minijets. We 
will show how the associated production of M ~ 1 — 2 GeV/c 2 dileptons from minijets with 
Pt > 2 GeV/c is comparable to the first order DY results at BNL Relativistic Heavy Ion 
Collider (RHIC) energies but becomes dominant at CERN Large Hadron Collider (LHC) 
energies, even up to masses M ~ 5 — 10 GeV/c 2 . We will also study the effects of nuclear 
modifications of the parton distributions, especially nuclear shadowing, to the dilepton rates. 

The remainder of the paper is organized as follows. In the next Section, we will calculate 
the dilepton fragmentation functions in the framework of pQCD. The connection and the 
difference between real photon fragmentation functions are discussed. We will derive both 
the first order result and the one with QCD evolution, including corrections to all orders 
in pQCD. In Sec. Ill, the dilepton fragmentation functions are convoluted with hard and 
semihard parton scattering cross sections to calculate the minijet-associated dilepton pro- 
duction in heavy ion collisions at both RHIC and LHC energies. Finally, a summary with 
some discussions on the implications to the dilepton production from a QGP is given in Sec. 



IV. 



II. DILEPTON FRAGMENTATION FUNCTIONS 

In this Section, we review the dilepton fragmentation functions of quark and gluon jets. 
We will work in an axial gauge so that interference terms in the final state radiation disappear 
in the leading logarithm approximation ||16| . 



A. Lowest order in pQCD 

Let us define z as the fractional light-cone momentum and q 2 as the virtuality of an off- 
shell parton as illustrated in Fig. |T](a). The differential cross section for a quark produced 
in a hard process with momentum scale Q to emit a dilepton with invariant mass M is, 



a dzdz e dq 2 dM 2 i 2irq 2 q ^ qK j 2ttM 2 
where is the fractional charge of the quark q iy <jq is the total cross section of the hard 
process, and 

*W*) = - z [l + (l-z)% (2) 
P^ e+e -(z) = z 2 + (l-z) 2 } (3) 

are the splitting functions for q — > jq and 7 — > £ + £~ which are similar to those of q — > gq 
and g — ► qq, respectively, except for the color factors. Integrating over the virtuality q 2 of 
the intermediate quark and the fractional momentum zi of one of the leptons, one has the 
QED dilepton fragmentation function of a quark, 

da 



D% qi {z,M^Q 2 ) = j Q M Jq 2 [d Z ^ 



do dz dze dq 2 dM 2 



One can see that D^^ q (z,M 2 ,Q 2 ) is similar to a virtual photon fragmentation function, 
except for a factor due to the extra QED coupling and the integration over the relative phase 
space of the leptons, 



D$i hi (z,M\Q>) = SL^Lrito^tf^cf). (5) 

For a real photon fragmentation function, the lower limit for the integration over q 2 in 
Eq. |4] is in principle given by the quark mass. For massless quarks, the infrared divergence 
in the lowest order has to be regulated by some cutoff of the hadronic scale. In the absence 
of a large mass scale, the QCD corrections to real photon fragmentation function have also 
to be regulated by some cutoff. The physics below the cutoff becomes nonperturbative. 
One has to introduce some initial conditions for the real photon fragmentation functions, 
either given by experimental data or by some model-dependent assumptions. The problem 
of dilepton production is different because the fixed invariant mass M provides a natural 
cutoff below which kinematic restrictions will terminate the processes. The QCD processes 
above this cutoff are in principle calculable to all orders. 

Since gluons are not directly coupled to photons and leptons, the dilepton fragmentation 
function of a gluon in the lowest order is, 

D^ /g (z,M\Q 2 ) = 0. (6) 

For later convenience, we define 

ln(Q 2 /A 2 ) 



k(M, Q) = In 
and 



ln(M 2 /A 2 ) 



(7) 



Qf\z) = e*kl + {l-z)\ (8) 



so that we can rewrite Eq. |3] as 

<J*,^,« = (£)'^m£)(*- -!»!«>(*). (9) 



B. First order contributions 



The first order contribution in pQCD to the dilepton fragmentation function of a quark 
comes from a gluon bremsstrahlung before the virtual photon production as shown in 
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Fig. 0(b). Remember now that z = Z\z^ is the fraction of the momentum, carried by 
the dilepton, of the initial quark before the gluon radiation. Defining the convolution of two 
functions as 



A <g> B(z) = f 1 ^A(z 1 )B(z/z 1 ), 

J z Z\ 



(10) 



it is straightforward to write down the first order dilepton fragmentation function, 



a 



2nql J z Z\ 
2 2 . ,M\ 2 



2W SW^W 8 "- 1 -"^™ ™ 



where, 



(11) 



Air 



Po = ll-2n//3, 



(12) 



A,hi(g 2 /A 2 ) 

is the running strong coupling constant with n/ quark flavors. The splitting function for 
q — > g<? in QCD is 



Pq^qg{ Z ) 



1 + z 2 
1 - z 



(13) 



The "+ function" here is introduced to include the virtual corrections to cancel the singu- 



larity from the soft gluon emission and to guarantee momentum conservation 1 16,17]. Other 



splitting functions we will use in the following are the standard ones [17] , 



Pq->gq( z ) 



Pg^qq( z ) 



Pg^gg( z ) 



Pq^qgi^ z n 



l-z) 



1 — z 
+ + z{\ - z) 

z 



.11 2n } 
12 ~ ~36~ 



)8(l-z) 



The convolution in Eq. [TT| can be easily done and it gives, 



(14) 
(15) 

(16) 



Qr\ z ) 



2 



_L {2[1 + (1 - zf\ ln(l - z ) + (2- z)z In z + z{2 - -z) 



(17) 
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To the first order in pQCD, the dilepton fragmentation function of a gluon is not zero 
anymore. From the diagram in Fig. |l](c), we have 

D£l /g (z, M\ Q 2 ) = £ /J dqf^L [ ^ W*0^/J^ M\ Q 2 ) 

\ 2 O 71 /f 2 o ^ n / 

" ~ H^) w {^-i-^)Y. p ^m®Qf\^ ( 18 ) 



,2W 3M 2 v A 2y /3 — x 
and the convolution in z can also be calculated explicitly in this case, giving 



2n f 

Mi 



1=1 

2n f 



= E e *Yz [ l {1 " z3) + z{1 ~ z) + 2(1 + z)z ln (19) 

As we will see below numerically, radiative corrections to any order will soften the QED 
fragmentation function of quarks and increase the fragmentation function of gluons. Because 
of the leading logarithm behavior of the radiations, the dependence of QCD corrections on 
the strong coupling constant is cancelled out so that they might become important to all 



orders. From Eqs. |9] and [TT|, we can see that the relative importance of the first order QCD 
correction to the QED fragmentation function is controlled by a Q-dependent factor, 

c ~ l ~^rr (20) 

where k is defined in Eq. [7[ For values of Q 2 not too large relative to M 2 , k is very small so 
that higher order corrections can be neglected. Only for extremely large values of Q 2 and 
thus k, C becomes comparable to 1. Then one has to include corrections to all orders. For 
our consideration here, Q 2 is in the order of p\ of the minijets. Thus, as we will show in the 
next Section, for most of the minijet production with px ~ 2 GeV/c, first order calculation 
of the dilepton fragmentation functions should be sufficient. 



C. Full QCD evolution 



Following the same steps as we have calculated the first order corrections to the dilepton 
fragmentation functions, we can calculate the higher order contributions. Here we neglect 

7 



the further splitting of the radiated soft gluons and quarks, and only consider those diagrams 
with a simple ladder structure in leading logarithm approximation. Therefore, the radiated 
soft gluons and quarks are always on mass-shell. The general form of the contributions with 
n radiations before the dilepton production can be derived as 

2 2 , (M 2 \ 



a 



2tt/ 3M 2 

'_2_y 



In 



A 2 



(21) 



D& /g (z,M*,C? 



a 



2 2 



2tt/ 3M 2 



In 



M 2 
A 2 



Jo, 



e K - (1 + k + 



K 



G {n \z) 



(22) 



where Qj n) (z) and G (n \z) can be calculated iteratively from the lower order results via 



2n f 



>!->qg ® Qi (z) + Pq^gq 



G (n_1) (; 



(23) 
(24) 



Since we know Q\ (z) (see Eq. |) and G {0 \z) = 0, we can in principle perform the above 
convolutions up to any order as we did for Q\ 1] (z) (Eq. 0) and G^(z) (Eq. |l|), and obtain 
the full QCD dilepton fragmentation functions as, 



D» h/qi (z, M\ Q 2 ) = £ i?W ( Zj m 2 , Q 2 ) 



n=0 

oc 



D DL/s (*, M 2 , Q 2 ) = £ DW /g (z, M 2 , Q 2 ). 



(25) 
(26) 



n=l 



For large values of n, evaluating the integration in the convolution analytically is obvi- 
ously too cumbersome. One method to evaluate the full fragmentation functions is to solve 
numerically a set of coupled evolution equations. 

Taking derivatives of D DL / q . g (z,M 2 ,Q 2 ) with respect to Q 2 and using 

dn 1 



(27) 



dQ 2 ln(Q 2 /A 2 )' 

and the definition of a s (Q 2 ) in Eq. |T^, we can derive from Eqs. pi|-|2"^| the following coupled 
evolution equations, 



S 



+ ^7^ P ™ ® Adl/^, M 2 , Q 2 ), (28) 
^(z, M\ Q 2 ) = £ P 9 ^ ® D DL/9i (,, M 2 , Q 2 ) 

+ ^2^ P ™ ® Pdl/ ^' Q2) - (29) 

These evolution equations are very similar to those of real photon fragmentation functions 
fnyiBfl and the parton distribution functions in a photon [Tj|. The only difference is that 
dilepton (or virtual photon) fragmentation functions with a given mass M have a definite 
initial condition, 



D DL/qi , g (z,M 2 ,Q 2 ) Q2=M2 = 0, (30) 

together with the boundary condition, 

D DL/g ^(z,M 2 ,Q 2 )\ z=i = 0. (31) 

The boundary condition simply means that the probability for the dilepton to take the whole 
fraction of momentum of the initial quark or gluon is zero after QCD evolution is taken into 
account. Since there is always a finite contribution to dD^/qJ dQ 2 from the QED term in 
the evolution equation Eq. |28], one can verify that Dm,/ qi (z, M 2 ,Q 2 ) must approach zero 
as l/ln(l — z) at z = 1 in order to satisfy the boundary condition at all Q 2 . For gluons, 
Dr>L/g{z, M 2 , Q 2 ) must go to zero faster than 1/ ln(l — z). 

The above evolution equations with the initial and boundary conditions can be solved 
numerically. The scale M 2 now sets the starting point of the evolution. We show the QCD- 
evolved dilepton fragmentation functions zDm,/^ g (z, M 2 ,Q 2 ) scaled by a common factor 
(a/27r) 2 (2/3M 2 ) ln(Q 2 /M 2 ) in Fig. | for M = 1 GeV/c 2 and Q = 5 GeV. Together, we also 
show the analytical results to the lowest and first order. It is clear that both the first order 
corrections and the full QCD evolution soften the fragmentation functions. The overall QCD 
corrections to the QED (or lowest order in pQCD) result are about 10%, except near z = 



and 1. Since a gluon does not have dilepton production to the zeroth order in pQCD, the 
dilepton fragmentation function of a gluon is one order of magnitude smaller than a quark. 
Because there are logarithmic divergences at z = and 1 for each order correction to the 
dilepton fragmentation functions, as can be seen in Eqs. [L?| and [19], every order becomes 
important so that one has to sum them together to get the full QCD result. This is why the 
full QCD-evolved fragmentation functions in Fig. ||] differ considerably from the first order 
results near z = and 1. To the first order, the fragmentation function of a quark is exactly 
proportional to the square of its fractional charge, ef. This charge scaling is only slightly 
violated at small z for large Q due to the gluonic contribution to the QCD evolution as seen 
in Eq. 



III. ASSOCIATED DILEPTON PRODUCTION 
A. Kinematical limits 

In this paper, we are interested in the dilepton production cross section integrated over 
the transverse momentum. Hence, we need only the dilepton fragmentation functions in- 
tegrated over z. As we have seen in the previous Section, the fragmentation functions 
diverge at z — 0. One must therefore introduce an infrared cutoff. Fortunately, for dilepton 
production, the invariant mass M provides a natural cutoff. 

Assuming that Q 2 and q 2 are the virtualities of the parton before and after the emission 
of a virtual photon with fractional momentum z, one can verify that the relative transverse 
momentum of the dilepton with respect to the original parton is 



k\. = z(l — z) 



M 2 q 2 



Q 



(32) 



z 1 — z 

Neglecting q 2 and requiring k\ > 0, we can see that M 2 provides a natural kinematical 
cutoff for z, 

z>z = M 2 /Q 2 . (33) 
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In principle, one could take into account these kinematical limits at every step of the ra- 



diation processes, as done in Monte Carlo approaches [|19|-|22|j. Although not shown here, 
this can be done analytically for the first order calculation of the dilepton fragmentation 
functions. One could also use the relative transverse momentum fcf, as the argument in the 
running strong coupling constant. This is, however, beyond the scope of our simple leading 
logarithm estimates in this paper. 

With the kinematical cutoff in Eq. |33], we can obtain the integrated dilepton fragmenta- 
tion functions, D^i qug {M 2 , Q 2 ), the probabilities for a quark or gluon to produce a dilepton 
with mass M within the interval dM 2 . The lowest and first order fragmentation functions 



can be obtained analytically by integrating Eqs. || [IT] and |1^ over z, 

a\ 2 2 , ,M 2 



< s ,(M*,Q=) = e?(^) 3^1n(^)(e«-l 

, Q 2 3 M 2 lto 
-M 2 >~2 + ~Q 2 ~2~Q 4 



Q 2 3 M 2 1 M 4 1 
21n &)-7T + 27^-o7^ , (34) 



rV\ _ 2/ «\ 2 2 . ,M\ 2 , „ . ,4 



M 2 2 


(e K - 




M 2 \ 


-(- 


Q 2 ) 


•-¥> 


♦(■ 



(35) 



,M\ 2tt 2 / Af 2 \ / M 2 \ , , M 2 , 

'4 m 2 m 4 \ ,q 2 

1 / M 2 \ / M 2 M 4 \l , . 

where the function g%{x) is defined as 

00 n 

toW = E (37) 



n=l ^ 



We plot in Fig. | the full QCD-evolved results L> DL/ g ii9 (M 2 , Q 2 ) as functions of M 2 at 
fixed Q= 4 GeV. The first order results D^l /q .(M 2 , Q 2 )+D^ L/qi (M 2 , Q 2 ) and Df^M 2 , Q 2 ) 
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are very close to the full QCD-evolved fragmentation functions with only a few percent 
difference through the whole M 2 range. As we have seen in Fig. |2], the full QCD-evolved 
fragmentation functions are enhanced at small z while depleted at large z as compared to 
the lowest order calculations. For small values of M 2 /Q 2 , QCD evolution is important, 
but the lower limit Zq of the z-integration is also small. Thus, the integrated full QCD 
fragmentation functions are almost the same as the first order results. At large values of 
M 2 /Q 2 , the lower limit zq is large, but the QCD corrections in any order are increasingly 
smaller. Therefore, in the whole range of M 2 , the first order calculation of the z-integrated 
dilepton fragmentation functions is a very good approximation. 



B. Dilepton production associated with minijets 

In the following, we consider dilepton production associated with minijets. In particular, 
we are interested in the differential rates of dileptons with rapidity Y = as functions of 
the invariant mass M. As a first approximation, we can assume the dilepton to be produced 
collinearly with the parent quark or gluon. Then the differential cross section can be writ- 
ten down in a straightforward manner by folding the z-integrated dilepton fragmentation 
functions D DL / q .(M 2 ,Q 2 ) and D DL / g (M 2 ,Q 2 ) with the 2 — ^ 2 subprocesses of minijet pro- 
duction. We can also neglect the contribution from the initial state dilepton radiation, since 



the rapidities of these pairs are typically large like those of the initially radiated partons p3 
One has to take into account that the dilepton pair can be produced by either one of the 
final state partons, and connect this to the correct normalization of the integrated minijet 
cross section Oj et . In this way, the basic formula for the associated production of dileptons 
with Y ~ from minijets in a AA collision at impact parameter b can be written as follows: 

dN^(b) 1 r» f/* 



dM 2 dY 
da 



1 r f s l A 

^ / d 2 r± dp 2 T dy x dy 2 X ifa/A(xi, Q 2 , r ± )x 2 f b/A (x 2 , Q 2 , |b - rj) 



Pq abed 
9,9,9 

*cd 



(s, t, u) [D DL/C (M 2 , Q 2 m J5(Y - yi ) + D DL/d (M 2 , Q 2 m JS(Y - y 2 ) J, (38) 
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where the produced (anti)quarks and gluons (i.e. minijets) have transverse momentum 
Po < Pt < and the kinematical range of rapidities 

- In ^5 - e^ 1 j < y 2 < In _ e ^ . (40) 
The momentum fractions of the initial state partons are denoted by 

^^(e^+e^ 2 ), (41) 
and the Mandelstam variables in the parton-parton level for the massless partons by 

s = x x x 2 s = 2p 2 T (l + cosh(yi - y 2 )), (42) 

* = -p 2 T {l + e m ~ yi ), (43) 
u = -p 2 T {l + e yi ' V2 ). (44) 

The cross sections da ab ^ cd /dt ~ 0(a 2 ) for the various partonic subprocesses can be found 



e.g. in Refs. |14| , f24f . The parton density of a nucleus by our definition is 

f a /A(x, Q 2 , rj_) = t A (r±)R a/A (x, Q 2 , r±)f a/N (x, Q 2 ), (45) 

where t A (r±) is the thickness function of the nucleus which is normalized to / d 2 r±t A (r±) = 
A. The parton distribution in a nucleon is f a /N(x,Q 2 ), and the ratio R a /A(x,Q 2 ,r±) for 
the nuclear modifications to the parton distributions is both scale and impact parameter 
dependent |Tl|p5| . In the following, we will approximate the impact parameter dependent 
ratio R a / A (x, Q 2 , r ±) by its averaged value, 

R a/A (x, Q 2 ) = jf d 2 r ± t A (r ± )R a/A (x, Q 2 , r ± ), (46) 

so that 

fa/ A (x,Q 2 ,r ± ) « t A ( ri _)tf m (x,Q 2 ), (47) 

where the effective parton distributions per nucleon in a nucleus is defined as 
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fa/{N)M 2 ) = Ra/A(x,Q 2 )f a /N(x,Q 2 ). (48) 



In this paper we will use the set 1 of the Duke-Owens parton distributions [£(J for 
f a /N(x,Q 2 ). We use the scale dependent nuclear modifications for R a /A(x,Q 2 ) as studied 
in [p7| . Especially, we assume that at the lowest scale Q = 2 GeV, gluons are shadowed 
by the same amount as the structure function in deeply inelastic £A scatterings. Note 
that the normalization of Eq. ^ can be checked by setting the M 2 -integrated fragmentation 
functions to unity and integrating over Y; this will give us 2crj e t(po, y/s) as expected when 
integrating over the inclusive 2 — > 2 scattering cross section. 

As usually in the case of pQCD calculations, there are uncertainties in choosing the 
momentum scales both in the parton distributions and the fragmentation functions. We will 
choose the scale entering the parton distributions to be the transverse momentum of the jets, 
Q = pt, for Duke-Owens parametrization set 1 with A = 0.2 GeV ||14|| . The scale Q max m the 



dilepton fragmentation functions represents the maximum virtuality of the final state parton 
before any radiation. As we have emphasized in this paper, the dilepton fragmentation 
functions at large fixed mass do not have nonperturbative contributions. Therefore, unlike 
the scale in the parton distributions, Q max in large mass dilepton fragmentation functions 
is not correlated with the choice of A. Examining the matrix elements of a + b^a + b + 'y* 
processes, one can find out that the scale entering the leading logarithm term is one of the 
Mandelstam variables, s, —t, —it, depending on the channel of the specific process. However, 
in Eq. [J^, we convolute the fragmentation functions with jet cross sections which include 
different channels and their interference terms. Therefore, Q max in Eq. is only an effective 
momentum scale. From Eqs. |42|444|, we know at least that Qmax > Pt- We will discuss the 



sensitivity to the choice of Q max when we present the results of our calculation. 

Note also that for the minijet production the lower limit po of the integration over pt is 
a parameter which determines the division between calculable "hard" and model-dependent 
"soft" processes. Most of the minijets are produced with p T ~ p ~ few GeV/c, and they 
are basically nonresolvable as distinct £r-clusters, even in hadronic collisions f28 |. The 
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phenomenological value of po depends on the model for cr so f t of soft processes, the par- 
ton distribution functions and the corresponding scale choice. Since these issues can not 
be addressed within pQCD, the possible range of values of po has to be determined phe- 
nomenologically, in connection with a model for the soft contribution cr so ft to the particle 
production in pp and pp collisions [O,O|2^J50[ . We will use here p = 2 GeV/c, as suggested 



and studied in detail in Ref. plj . Although already exactly calculated for inclusive jet pro- 



duction fl32| , the 0(a z s ) contributions to the lowest order parton cross sections are simulated 



here by an overall factor K ~ 2. Clearly, the parameter po depends also on the size of the 
next-to- leading order terms. We want to point out, however, that the cross section for the 
associated dilepton production in Eq. [38] is much less sensitive to the choice of po than the 
minijet cross section itself. For Q max = Pt, the dilepton fragmentation functions vanish for 
M > pt- Whenever M > p , M takes over as an effective cutoff in the integration over p? in 
Eq. |38|. Therefore, the cross section for the associated dilepton production does not depend 
on the exact choice of po at large M. 

The symmetrized formula of Eq. [SB] can be somewhat simplified by considering all the 
possible pairs of partons in the initial and final states, (ab), (cd). By changing the integration 
variables yx,2 into — 1/2,1 appropriately in the other half of the expression, and by using the 
t, M-symmetries of the subprocess cross sections, a Y «-> —Y symmetric formula can be 
written down. Especially, at Y = we get: 

dN%A /ict (b)\ _ ... /-/4 



dM 2 dY 



Y- 



rs/'i I I 

= 2T AA (b) dp 2 T dy 2 J2 ; , x ; , x ^ifa/(N)( x u Q 2 )x2fy {N} (x 2 , Q 2 

J Po /ab) 1 1 °ab J- -T Ocd 



(ab) 
(cd) 



{D DL/C (M 2 , QlJ-^ (s, t, u) + D DL/d (M 2 , QLx)^ (a, u, t) 



(49) 

2/1=0 



where Taa{^) = J d 2 r±tA(r±)tA.(\h — r±\) is the nuclear overlap function of the two colliding 
nuclei. 

As discussed in the previous Section, the first order results in Eqs. |3^ - p6| are a good ap- 



proximation for the full z-integrated dilepton fragmentation functions, which is the approx- 
imation we shall adopt in what follows. The results from Eq. (i^ with nuclear modifications 
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to the parton distributions are shown in Fig. f| (solid curves) for = 200 AGeV and 6400 
AGeV, respectively. In the figure, we have compared the minijet associated production of 
dileptons to the lowest order differential cross section of the direct Drell-Yan process (dashed 



curves 



dNZ(b) 



dM 2 dY 



= TAA(b)^J2 e 2 qi Lfi m (*i, M 2 )x 2 f£ /{N) (z 2 , M 2 ) 

i 

{x u M 2 )x 2 f* m {x 2 ,M 2 )l (50) 



y=o 



where Xi j2 = M/\fs at Y — 0. We have chosen the scale in the parton distributions as 
Q = M. To simulate the first order pQCD contributions to the DY cross section |||7|], we 
multiply Eq. |50] by an overall factor Kdy ~ 2. Note that since we have used the Duke- 
Owens parton distributions, which extend only down to Q = 2 GeV, the results for direct 
Drell-Yan cannot really be trusted much below M = 2 GeV/c 2 . 

To study the sensitivity of the minijet associated dilepton production to the choice of 
the scale Q max i n the fragmentation functions, we plot in Fig. |] the results of Eq. ^ for 
both Q ma x = Pt and 2pT- It is apparent that the results are relatively sensitive to the 



choice of Q max . As we can understand from Eqs. p4] - |36| , the difference between the two solid 
curves is due to the fact that the ^-integrated fragmentation functions are proportional to 
ln 2 (Q max /M 2 ). Due to the kinematical restriction M < Qmax, changing Q max = 2px to px 
also effectively doubles the lower limit of the integration over p^ for fixed M in Eq. |4£]. This 
is the reason why the two solid curves have different slopes. As one of the main purposes 
of this paper, Fig. |] demonstrates how the relative contribution of the dileptons associated 
with minijets in the range 1 <M< 10 GeV/c 2 changes with increasing energy as compared 
to the direct Drell-Yan production. Even after taking into account the uncertainties due 
to different choices of Q max , it can be seen clearly that at RHIC energy, yfs = 200 AGeV, 
dileptons from the bremsstrahlung of minijets are comparable to the direct Drell-Yan at 
M<2 GeV/c 2 . However, when going up to TeV energy range, dileptons associated with 
minijets become more important, and dominate the Drell-Yan at LHC energy, y^i = 6400 
AGeV, even up to masses M ~ 10 GeV/c 2 . Qualitatively, our results are similar to the 
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minijet-associated photon production in Ref. where real photon fragmentation functions 
in the lowest order are considered. 

To demonstrate the effects of parton shadowing and antishadowing, we plot in Fig. |5| 
the results calculated with (solid) and without (dashed curves) nuclear modifications of the 
parton distribution functions. We can see that nuclear shadowing depletes the Drell-Yan 
dileptons relatively more than the dileptons from the minijets. The basic reason for this is 
that Drell-Yan pair production dN^X m Eq. |5^ as a function of M = xi^y/s probes the 
(anti)quark distributions directly, at least in the lowest order. Furthermore, the antiquark 
shadowing does not vary strongly with the scale M, as has been experimentally measured 



34,35]. On the other hand, in the minijet-associated dilepton production, we have to in- 



tegrate the contribution over the whole range of x. In addition, we also have to integrate 



over the scale Q = pr- The gluon shadowing [27] we used here has stronger Q dependence 



than the (anti)quark. Therefore, the net effect of the nuclear modifications of the parton 
distributions to the minijet-associated dilepton production remains relatively small even at 
TeV energy range. 

IV. SUMMARY AND DISCUSSION 

In this paper, we have studied minijet-associated dilepton production in ultra-relativistic 
nuclear collisions. We calculated both the first order approximation and the full pQCD 
evolution of the dilepton fragmentation functions of produced partons. The dilepton pairs 
from the fragmentation of minijets are found to be comparable to direct Drell-Yan at RHIC 
energy for small invariant mass M ~ 1-2 GeV/c 2 . At LHC energy, the associated dilepton 
production becomes dominant over a relative large range of the invariant mass. These 
dileptons plus the direct Drell-Yan pairs would constitute part of the background to the 
dilepton production from a QGP and its pre-equilibrium stage. Other background includes 



dileptons from final hadronic rescatterings p6|j37|] and the decay of charmed hadrons [33,38 



It is also straightforward to calculate the pt distribution of the associated dilepton pairs in 
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our fragmentation function approach. Since one has to convolute the dilepton fragmentation 
functions in z together with the pr distributions of the jets, we expect the resultant px 
spectrum of these dileptons to be softer than the px spectrum of the jets. Therefore, the 
dileptons associated with minijets should have smaller pt relative to the direct Drell-Yan 
pairs which have a high p? tail like that of the produced jets. Since thermally produced 
dileptons in a QGP also have relatively small pt as compared to Drell-Yan Q, minijet- 
associated dileptons thus pose a more intangible background. 

In calculating the dilepton fragmentation functions, we have assumed leading logarithm 
approximation so that we can include contributions from all orders in pQCD. However, 
the higher order corrections are small and the first order results are sufficient enough for 
our estimates of the minijet-associated dilepton production. The largest uncertainty in our 
calculation is the choice of the momentum scale Q max used in the dilepton fragmentation 
functions. Since the correct scale in a matrix element calculation is channel-dependent, 
we used only an effective scale choice in the fragmentation functions to convolute with the 
minijet cross sections. We evaluated the dilepton spectrum for two choices of the scale, 
Qmax = Pt, %Pt- However, the results with Qmax = Pt should give us the lower bound of 
the associated dilepton production. Another notorious uncertainty of the px cutoff p in 
minijet-related problems is greatly reduced here due to the kinematic restriction M < Q max - 
For Qmax = Pt, the px cutoff is replaced by M whenever M is larger than p . 

The abundance of dileptons associated with minijet production at high energies is mainly 
due to the large gluon-related minijet cross sections and the high initial gluon densities inside 
the colliding nuclei. This should have important implications for the dilepton production 
in the pre-equilibrium stage of the quark gluon plasma. As pointed out recently |39|-|41|, 
the parton system is not at all in chemical equilibrium when initially produced in the ear- 
liest stage of high energy nucleus-nucleus collisions. Because of the small cross sections 
for (anti) quark production, the initial parton system is dominated by gluons and is quark 
deficient as compared to an equilibrated QGP. Studies P§|-fn|] also suggest that the parton 
system thus produced may not be able to achieve chemical equilibrium before hadronization. 
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In this case, dilepton production through qq annihilation should be severely suppressed. On 
the contrary, dilepton production from gluon fragmentation could become relatively im- 
portant for a gluon dominated system, since gluon-related cross sections of small angle 
scatterings are about 9/4 larger than the quark. Even though the dilepton fragmentation 
function of a gluon is about one order of magnitude smaller than a quark, a gluon density at 
least about 5 times higher than the quark could easily compensate the small fragmentation 
function and make the gluon associated dilepton production important. 
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FIGURES 

FIG. 1. Illustration of the diagrams of (a) the lowest order, (b) the first order contributions in 
pQCD to the dilepton fragmentation functions of quarks and (c) gluons. The dashed lines present 
the associated hard processes with momentum scale Q. 

FIG. 2. The QCD-evolved (solid), the lowest order (dot-dashed) and the first order (dashed) 
approximations of dilepton fragmentation functions zD DL/ / a (z, M 2 , Q 2 ) of a u-quark and a gluon, 
for M = 1 GeV/c 2 and Q = 5 GeV. A factor (a/2vr) 2 (2/3M 2 ) ln(Q 2 /M 2 ) is divided out. 

FIG. 3. The z- integrated dilepton fragmentation functions DrjL/ a (M 2 , Q 2 ) for a u-quark (solid) 
and a gluon (dashed) as functions of M 2 at fixed Q = 4 GeV. A factor (a/27r) 2 2/3M 2 is divided 
out. The curve for gluon fragmentation is multiplied by 10. 

FIG. 4. Mass spectra of minijet- associated (solid curves) and Drell-Yan (dashed) dileptons at 
Y = in central Au + Au collisions at yfs = 200 and 6400 AGeV. The two solid curves correspond 
to two choices of the scale Q max = Pt and 2px in the dilepton fragmentation functions. Parton 
shadowing is included in the calculations. 

FIG. 5. Mass spectra of the minijet-associated and Drell-Yan dileptons in central Au + Au 
collisions at y/s = 200 and 6400 AGeV, with (solid) and without (dashed) parton shadowing. 
For the associated production, the scale in the dilepton fragmentation functions is chosen to be 
Qmax = 2pr- 
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